Fidelity susceptibility, scaling, and universality in quantum critical phenomena 
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We study fidelity susceptibility in one-dimensional asymmetric Hubbard model, and show that the 
fidelity susceptibility can be used to identify the universality class of the quantum phase transitions 
in this model. The critical exponents are found to be and 2 for cases of half-filling and away from 
half-filling respectively. 
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Quantum phase transitions (QPTs) at zero tempera- 
ture are characterized by the significant change in the 
ground state of a many-body system as a parameter 
A in the system Hamiltonian H{X) is varied across a 
point Ac This primary observation enlightens peo- 
ple to explore the role of fidelity, a concept emerging 
\ from quantum information theory [3| , in the critical phe- 
nomena [1, 3. Since fidelity is a measure of similarity 
between states, a dramatic change in the structure of the 
ground state around a quantum critical point should re- 
sult in a great difference between the two ground states 
on the both sides of the critical point. Such a fascinating 
prospect was firstly confirmed in the ID XY model where 
the fidelity shows a narrow trough at the phase transition 
point [1, 01 • From then on, the fidelity was further used 
to characterize the QPTs in fermionic Q and bosonic 
systems @. As fidelity is purely a quantum information 
concept, an obvious advantage is that it can be a promis- 
ing candidate to characterize the QPT 0, i, S [13, [H El 
because no a priori knowledge of the order parameter and 
the symmetry of the system is needed. Therefore, these 
works established another connection between quantum 
information theory and condensed matter physics, in ad- 
dition to the recent studies on entanglement in QPTs 




The fidelity actually reflects the response of the ground 
state to a small change of the driving parameter. Za- 
nardi et al. introduced the Riemannian metric tensor Q 
inherited from the parameter space to denote the leading 
term in the fidelity, and argued that the singularity of 
this metric is in correspondence with the QPTs. While 
You et al introduced another concept, so-called fidelity 
susceptibility (FS) [£], and established a general relation 
between the leading term in the fidelity and the struc- 
ture factor of the driving term in the Hamiltonian. This 
relation implies that the fidelity may not have singular 
behavior in those transitions of infinite order, such as 
Kosterlitz-Thouless (KT) transitions [19]. 

In this work, we study the FS in ID asymmetric Hub- 
bard model (AHM) [2l|, and show that the FS can be 



'Electronic address: ' sjgu@phy.cuhk.edu. hk| 



used to characterize the universality class [205 in quantum 
critical phenomena. The intrinsic relation between the 
FS and the Landau's symmetry-breaking theory (LSBT) 
is firstly clarified by a simple QPT occurred in a well- 
studied ID transverse-field Ising model. Then we mainly 
focus on the critical behavior of the FS in the ID AHM. 
Since the AHM can be used to describe a mixture of two 
species of fermionic atoms in an optical lattice, which 
is able to be realized by recent experiments on the cold 
ator ns |22i. the model itself is of current research interest 
[H Him, m, m m. we find that the critical expo- 
nents of the FS take the value of and 2 for cases of 
half-filling (n ~ 1) and away from half-filling (n = 2/3) 
respectively. 

To begin with, we consider a general Hamiltonian of 
quantum many-body systems, i.e. 

H{\ h) = Ho + XHi + hM, (1) 

where Hj is the driving Hamiltonian with the strength 
A, and M is a potential order parameter and h is the 
corresponding external field. Without loss of generality, 
we first set h = 0. Following Ref. 0], the fidelity is 
defined as the overlap between two ground states |4'o(A)) 
and |^'o(A -I- i5A)), that is 

F(A,<5A) = |(1'o(A)|vI/o(A + a))|. (2) 

Then the FS is just the most relevant term in the fidelity, 
and mathematically is related to the structure factor of 
the driving term Hi Q, which denotes the fluctuation 
caused by the driving parameter. For_ example, if we 
extend the fidelity to the thermal state [3], the FS is just 
the specific heat or the magnetic susceptibility [1, Q if we 
choose the driving parameter as temperature or magnetic 
field respectively. 

Compared with ordinary phase transitions, we can also 
extract two exponents (denoted as 01,7) from the FS in 
quantum critical phenomena if we choose the driving pa- 
rameter as A and h (if the order parameter is known) 
respectively (Here, the only condition is that the QPT 
should belong to the type of Landau's transition, other- 
wise a = and 7 is not well defined). Then the fidelity 
susceptibilities driven by two terms in the Hamiltonian 
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respectively, around the critical point Ac in the thermo- 
dynamic limit. As a simple application, we take the well- 
studied model, i.e. ID transverse-field Ising model, as an 
example. 



where cr is Pauli matrix. In Ref. [J], it was obtained 
a = 1. On the other hand, if we consider h as the driv- 
ing parameter, we find 7 = 7/4, which is the same as 
the 7 of the 2D Ising model [23|. Take into account the 
exponent (3 — 1/8 of the order parameter cr^, we then 
have a + 2P + J = 3, which is slightly different from the 
usual a + 2(3 + J = 2 in 2D Ising model. We interpret 
this difference as one more differentiation is made in the 
FS of the ground state than the specific heat at finite 
temperatures. That is, the FS is related to the second 
order derivative of the ground state energy with respect 
to the driving parameter while for the specific heat, 
it is simply dE{T)/dT where E{T) is the internal energy. 
Therefore, the phase transition here still belongs to the 
same universality class of the 2D Ising model. As a brief 
conclusion, the relation between the FS and the LSBT 
is straightforward. Once the driving term and order pa- 
rameter are given, the universality classes are simply de- 
scribed by the critical exponents of FS. 

The LSBT is established on the order parameter, 
whose non-vanishing behavior results from the broken 
symmetry and long-range order. For KT transitions, 
both broken symmetry and long-range order are absent, 
hence no local order parameter is concerned. Consider 
again the original definition of the FS, i.e. 
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where |*„(A)} satisfies i7(A)|*„(A)) = £'„|*„(A)} and 
defines a set of orthogonal complete basis in the Hilbert 
space. For the KT transition, despite of the vanishing 
energy gap, there is still no singularity in the FS as ma- 
trix elements (^'„(A)|i?/|\['o(A)) also vanish at the same 
time. However, the appearance of the power-law decay 
behavior describes the stronger fiuctuation around the 
critical point. This point directly leads to that the FS, 
which also denotes the fluctuation of the driving term, 
might reach a maximum near the critical point, though 
the maximum point might not be the critical point, as 
has been observed in the ID Hubbard model 0. 

To confirm this physical intuition, we now focus on the 
ID AHM, whose Hamiltonian reads 
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FIG. 1: (color online) The schematic phase diagram of the 
AHM. LEFT: the phase diagram defined on the U — t plane 
which shows a KT transition at half-filling and Landau's tran- 
sition away from half-filling. RIGHT: the phase diagram de- 
fined on the n — t plane, in which the transition along the 
middle line is of KT type and the thin solid lines at both 
sides are of Landau's type. The phase diagram has a mir- 
ror symmetry about the line n = 1 due to the particle hole 
symmetry in the model. 



where ^ and Cj^o-, o' —1, 1 are creation and annihilation 
operators for electrons with spin a at site j respectively, 
Ua = cJ-Cct, to- is cr-dependent hoping integral, and U 
denotes the strength of on-site interaction. In this model, 
the Hamiltonian has U(1)(8)U(1) symmetry for general to-, 
and the atoms number 7V| = X]j '^i.T' -^i = '^i.i 
conserved respectively. The total number of atoms is 
given by N = + Ni, and the filling factor is n = N/L. 
For simplicity, we reset t = ti/t^, and U to be U/t^. 

The schematic phase diagram of the AHM is shown 
in Fig. [TJ which can be understood from its two lim- 
iting cases, the Hubbard model [30] (i| — ti) and the 
Fahcov-KimbaU (FK) model [HI, 13^ {ti =0). At half- 
filling, both the Hubbard model and the FK model are 
in a spin-density- wave state. The difference is that in the 
Hubbard region, the system renormalizes to the Heisen- 
berg fixed point, while in the FK region, it belongs to the 
Ising fixed point. The QPT occurred between these two 
classes belongs to the KT type in the ID system because 
the correlation function at both side is of power-law decay 
[2l| and no local order parameter is well defined. While 
away from half filling, the system becomes an ideal con- 
ductor and is in the state of density wave in the Hubbard 
region, but it is in a phase separation in the FK region. 
In the phase separation region, the translational sym- 
metry is broken, and the down-spin electrons congregate 
together; then (n^) plays a role of the order parameter. 
So the phase transition is of Landau's type t27i.j28j. 

In order to quantify the change of the ground state 
during the evolution of t, we define the fidelity as 
F{t,5t) = \{-9{t)\-^{t + 5t))\. The corresponding FS 
is XF{t){t) = — 2hm5f^olnF(<,(5t)/(5t^. To avoid the 
ground state level crossing, we choose the periodic or 
antiperiodic boundary conditions for systems with 4Z -I- 2 
or 41 electrons respectively. We first look at two spe- 
cial cases, i.e. the FS at half-filling (n — 1) and away 
from half-filling (n — 2/3), and both with a given inter- 
action U = 10. The numerical results of different sys- 
tem sizes are presented in Fig. [21 For both cases, the 
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FIG. 2: (color online) The scaling behavior of the FS as a 
function of t for the cases of n = 1 (LEFT) and n = 2/3 
(RIGHT). Here U= 10. 
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FIG. 3: (color online) The finite size scaling is performed 
for the case of power-law divergence for the case of U — 30, 
n = 2/3, and system sizes L = 6,9, 12, 15. The FS, consider 
as a function of the system size and the driving parameter is 
a function of L"{t — imax) only. Here the critical exponent is 
v ~ 2.65. The inset denotes the scaling behavior of XF{t ~ 
imax). The straight line is of slope 1 in logarithmic scale, and 
^ ~ 5.3. 



FS reaches a maximum point at a certain position tmax- 
The difference is that for n = 2/3 case, XF{t)it = ^max) 
diverges with increasing system size; while for rt = 1, 
XF{t)(t = ^max) oc L. The former behavior clearly de- 
notes a Landau's transition, and the latter is KT transi- 
tion. 

To study the critical behavior around the critical point, 
finite scaling analysis is needed. According to the defi- 
nition of the critical exponents in Eq. ([3]), we introduce 
the following scaling behavior for the FS 
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FIG. 4: (color online) The similar finite size scaling is per- 
formed for the case of KT transition occurred for various sys- 
tem size L = 6,8,10,12 at n = 1 and C/ = 30. The FS, 
consider as a function of the system size and the driving pa- 
rameter is a function of L" {t — fmax) only. Here the critical 
exponent is ~ —0.25. The inset denotes the scaling behavior 
of X^ii = imax). 



where A, B are constants independence of L and t. 
Such a finite size scaling leads to that the rescaled FS 
[XF(t){t = tmax) - XF(t){t)]/xFit){t) is a simple function 
of the rescaled driving parameter i''(i — <,nax)- This func- 
tion is universal and does not depend on system sizes, as 
shown in Fig. [3] for the case of J7 = 30, in which nu- 
merical results obtained from various system sizes fall 
onto a single line. On the other hand, if t = fmax which 
approaches to tc like tmax — tc oc L~^, the maximum 
value of the FS diverges with increasing system size as: 
XF(t)it = imax) C!C . Then the exponent ly together 
with fj, determines the critical exponent a in Eq. ([3]). 
For the present case, we find v ~ 2.65 and fi ~ 5.3, 
hence a = fi/v = 2. That is, around the critical point, 
the FS for the case of n = 2/3 scales like 



XF{t){t) 



1 



L 



(8) 



which clearly differs from the Ising model. For the Ising 
model [Eq. Q], only Z2 symmetry is broken when the 
phase transition occurs; while in the AHM, the transla- 
tional symmetry is broken in the phase separation region. 
So they belong to different universality classes. 

The KT transition occurs at half- filling n = 1. The 
corresponding finite scaling analysis for the case U — 30 
is presented in Fig. |4l The maximum point of the FS 
is proportional to the system length. This is consistent 
with our aforementioned understanding. On the other 
hand, the rescaled FSs for various system sizes fall onto 
a single line, which is a function of L'^(t — imax) with 
critical length exponent 1/ = —1/4. We find that the FS 
around the maximum point like: 
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around tmax- As expected, there is no singularity in 
XF{t)it)- Clearly, the maximum behavior becomes weak 
with the increasing system size, as we can also find in 
Fig. [21 Then the FS becomes flat in the FK region, the 
relative larger FS is due to the power-law behavior the 
correlation function. On the other hand, motivate d by 
the KT transition occurred in quantum XY model [33|, 
we infer that the very steep decreasing point of the FS 
in Fig. [2] is more close to the critical point. Therefore, 
we perform 1/L finite scaling analysis for the minimum 
point of dxF{t) {t) /dt, and find ~ 0.308, 0.313,0.317 for 
U = 10, 20, 30 respectively. The results are very close to 
those obtained by density matrix renormalization group 
method 

A similar analysis can be carried on for other filling 
conditions. The power-law divergence of the FS always 
exists in other filling conditions except when n = 1. On 
the other hand, due to the particle-hole symmetry in the 
AHM, the FS takes the same value for n and 2—n (n < 2), 
which satisfies the same scaling behavior. Therefore, the 
phase diagram in the left picture of Fig. [T] has a mirror 
symmetry about the line n = 1. Take into account the 
fact that even a single hole doping might lead to the 
instability of the density wave state in the infinite U limit 
[131, the KT transition only happens at the half- filling 



condition. So the circle point in Fig. [T]is expected to be 
a quar-critical point in the phase diagram, and the FS 
just signals the transition type along the critical lines. 

In conclusion, we have shown that the FS, as the lead- 
ing term in the fidelity between two ground states sep- 
arated by a slightly difference in parameter space, can 
be used to characterize the universality class in quan- 
tum critical phenomena. Since the FS is related to the 
structure factor of the driving term in the Hamiltonian, 
the relation between the LSBT and FS is linked up. Its 
critical exponent then is naturally suitable for the clas- 
sification of universality. We elucidate this relation by 
the simple QPT occurred in the ID transverse-field Ising 
model. Furthermore, despite no singularity appearing in 
the FS when crossing a KT transition point, the stronger 
fluctuation might makes the FS reach a maximum near 
to the critical point. We then studied the FS in the ID 
AHM, and shown that the FS can help us to identify 
both types of phase transition in this model. The critical 
exponent a for the Landau's type transition is calculated 
with finite size analysis, and is found to be a = 2 for 
n = 2/3 case. While for the KT transition, a = 0. 

Note added. Recently, the work on the scaling behavior 
of the FS in other model appeared [ll|- 
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